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ABSTRACT 


This thesis reviews the technique established to clear channels in 
the Power Spectral Estimate by applying linear combinations of well 
known window functions to the autocorrelation function. The need for 
windowing the autocorrelation function is due to the fact that the true 
autocorrelation is not generally used to obtain the Power Spectral 
Estimate. When applied, the windows serve to reduce the effect that 
modifies the autocorrelation by truncating the data and possibly the 
autocorrelation has on the Power Spectral Estimate. 

It has been shown in previous work that a single channel has been 
cleared, allowing for the detection of a small peak in the presence of a 
large peak in the Power Spectral Estimate. 

The utility of this method is dependent on the robustness of it on 
different input situations. We extend the analysis in this paper, to 
include clearing up to three channels. We examine the relative 
positions of the spikes to each other and also the effect of taking 
different percentages of lags of the autocorrelation in the Power 

Spectral Estimate. 

This method could have application wherever the Power Spectrum is 
used. An example of this is beam forming for source location, where a 
small target can be located next to a large target. Other possibilities 
extend into seismic data processing. As the method becomes more auto- 
mated other applications may present themselves. 
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CHAPTER 1 


INTRODUCTION 

Autocorrelation estimates are adversely effected when truncated 
data sets are used (Blackman and Tukey, 1958). To overcome the bias due 
to the data truncation, the autocorrelations are themselves often trun- 
cated, producing sidelobes around peaks and the corresponding power 
spectral estimates. Figures 1.1 thru 1.3 demonstrate this effect. 
Figure 1.1 is a power spectrum with a single spike residing in channel 
31. When inverse Fourier transformed this spike produces a function 
domain time series shown in Figure 1.2. This function domain time 
series is modified by taking only a limited number of samples from the 
original series corresponding to the multiplication by a rectangular 
window. This modified or truncated time series is Fourier transformed 
so that we witness the effect of this modification (Figure 1.3). The 
spike remains prevalent in channel 31 but now the "ringing" sidelobes 
pervade throughout the spectrum. These side lobes can override and mask 
smaller peaks in the vicinity of the dominant peak. The application of 
windows also has the general effect of smoothing the power spectral 
estimate which may also have a deleterious effect on a small peak in the 
vicinity of a larger peak in the spectrum. 

There has been much work done with the utilization of window 
functions to reduce the effects of truncation. Harris (1978) details 
the use of windows in harmonic analysis, and describes a number of 
variables with which to measure the performance of each. Specific 
application of these windows has been used to modify an autocorrelation 
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function prior to power spectral estimation. This technique is de- 
scribed by Blackman and Tukey (1958) and Robinson (1980). The 
autocorrelation being windowed is modified such that when transformed 
the effects of this modification are seen in the power spectrum. These 
effects are typically seen as reduced sidelobes associated with dominant 
peaks in the power spectrum. Still there remains sidelobes, although 
they are much smaller than when no window used. 

Because of these effects, Smith (1985) developed a technique in 
which a linear combination of window functions were implemented to clear 
a single channel in the power spectral estimate. Although the technique 
was established, its implementation was limited to two input data 
situations, which will be reviewed. As a consequence, this thesis 
expands on the work done by Smith by examining the performance of this 
method on different input data sets. A review of the mathematical 
background and theory of the technique is given in order to establish a 
foundation for this work. To understand the procedures involved in 
utilizing this method, it is beneficial to understand the theory it is 
based on. In light of this, a review of the techniques and algorithms 
used in performing these calculations is presented. Finally, the input 
data and results of the calculations performed on these data are pre- 
sented and analyzed. 

This method has potential applicability in a number of fields which 
involve power spectral estimation. Among these is beam forming, where 
the objective is to locate a small target in the presence of a larger 
target. In seismic data processing we will, in future research, inves- 
tigate the influence of a dominant signal in the power spectrum, 
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specifically 60 Hz noise. The purpose of this is to determine by 
channel clearing whether we can detect reflection signals near a domi- 
nant event such as 60 Hz noise. 

The mathematical notations used in this paper will be similar to 
those Smith (1985) used in his work. The symbol w with or without 
subscripts will refer to a weighting function or a window in the func- 
tion domain (t) . The Fourier transform of w(t) is represented by W(s) 
where (s) is the transform domain variable. Likewise, the function 
domain autocorrelation is represented by g(t) and its Fourier transform 
by G(s) . 
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CHAPTER 2 


TECHNICAL BACKGROUND 


In order to understand the need for channel clearing in the Power 
Spectral Estimate, we need to review the concept and purpose of the 

autocorrelation function and windows. Given a discrete input time 
series with M terms 


(x(0) ,x(l) ,x(2) x(m-l) 

where 

x(j)-0 for j<0 and j>M-l, 


the autocorrelation is defined as follows: 

MH-| 

g(J T )- 2] x( k +j) x (lc) (1) 

k=o 

where j is a lag index and T is the sample increment. The 
autocorrelation function, as described by equation (1), is a shift, 
multiply and add operation which has the property of being non-negative 

definite (Robinson 1980). It is also commutative, which is expressed by 
the following equality: 

j —| M-j— 1 

^Tx(k+j)x(k)- ^ x(k)x(k- j ) 

k=° k=o 


( 2 ) 


5 


We can show this easily by expanding the summation operation as 
follows. We write the input series 

x(0) x ( 1 ) x(2) x(3) x(m-l) 

x(0) x(l) x(2) x(3) x(m-l) (3) 

where in this case j-0, meaning there is no time shift or we are 
calculating the zero lag value of the autocorrelation. The result of 
applying (2) to (3) gives 

g(0)-x 2 (0)+x 2 (l)+x 2 (2)+ 2 (3) + x 2 (m-l) . (4) 


We recognize equation (4) as the total energy of the series given 
in (3), because by definition (Robinson (1980) the total energy of a 
signal is given by 


m -1 

E- 2 |xf\ (5) 

n=o 


The lag value will now be shifted by 2 units to the right. 
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x(0) x(l) x(2) x(3) x(m-l) (6a) 


x(0) x(l) x(m-3) x(m-2) x(m-l) (6b) 


or 


g(2)-(x(0)x(2))+(x(l)x(3))+. . . . (x(m-l) ) . (7) 


Similarly a shift of 2 units to the left yields 


x(0) x(l) x(m-3) x(m-2) x(m-2) (8a) 


x(0) x(l) x(2) x(3) x(m-l) (8b) 


therefore 


g(-2)-(x(2)x(0))+(x(3)x(l))+ (x(m-l)x(m- 3) ) (9) 


Now we can see from equations (7) and (9) that 
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g(-2)-g(2) 


and in general , 


g(j)-gt-j) 


( 10 ) 


showing the autocorrelation function is symmetric about the zero lag 

U-0). 

It was shown earlier that the zero lag value of the autocorrelation 
of a signal was equivalent to the total energy of the signal. Another 
quantity calculated from the autocorrelation is based on the interrela- 
tion of the autocorrelation function and the power spectral estimate and 
the effects windowing has on each. Bracewell (1978) , through the use of 
the Autocorrelation Theorem and associated derivation, defines the 
relationship between the autocorrelation function and the power spec- 
trum. These two quantities are Fourier transform pairs where the 
Fourier transform of a function f(x) is defined as 


r °° 

F(s)— J f (x)e-»-2tf’5* dx. (11) 


Similarly, the Fourier transform of F(s) is defined as 
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f(x) 




F(s)e 


\lirsx 


ds. 


( 12 ) 


The quantity given in equation 11 is the result of what Harris 
(1978) refers to as a decomposition of a signal over a basis set. The 
Fourier transform specifically decomposes the signal set (f(x)) into the 
sums of sines and cosines. This complex arithmetic transform is defined 
for the discrete functions by the following: 


F(s)- T/T(0) 2 f(x)cos(2 7Tsx/M) - 


(13) 


iT/T(0) ^ f (x)sin(2 irsx/M) 


where 


M - the total number of samples 


T - the sample interval 


s - a multiple of the harmonic frequency counter, a multiple of 
1/T(0) 


T(0) - the length of the data set 
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x - the sample counter 


f(x) — input data samples 


Using the Euler relation 


iO 

e 


— cos# -isin# 


( 14 ) 


we can rewrite equation 13 as 


F(s)- 2 f (x)exp( -2rrsx/M) . 


(15) 


A discussion of the applicability of equations 11 and 15 will be 
given later in this chapter. The sums in equation 13 are over the 
aforementioned basis set. Through this decomposition or Fourier trans- 
form we obtain a representation of the signal in terms of frequency 
based on the component of each monotonic cosine and sine. The f(x) and 
F(s) in equations 11 and 12 are referred to as Fourier transform pairs 
and are represented by the following notation 
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f (x) 3) F(s) . 


Therefore , in agreement with the notation given in the introduc- 
tion, the Fourier pair, autocorrelation and power spectrum, are repre- 
sented by the following. 


g(x) G(s) . 


The calculation of the Fourier transform as defined in equation 11 
is for a continuous function of x, evaluated form minus infinity to plus 
infinity. Due to the reality of not being able to measure an infinite 
signal, we are confined to analyzing a finite data set. Also, we are 
not interested here in evaluating a continuous function but rather, a 
sampled function. The above discussion suggests the relationship 
between equations 11 and 15, i.e., equation 15 is the finite sampled 
function domain equivalent of equation 11. Brigham (1974) extensively 
discusses sampling theory as related to the Fourier transform and also 
discusses the algorithms developed for the Fast Fourier Transform, which 
is implemented in the work presented in this paper. Inherent in the 
transform of this sampled- finite data set, are artifacts of the numeric 
process. These artifacts of the transform process can cause the 
deleterious consequences shown in Figure 1.3. Specific to this, Harris 
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(1978) discusses the concept of spectral leakage. When sampling a 
finite data set with T sample interval, we obtain a total length nT, 
where n is the number of samples. If there are frequencies contained in 
the signal not contained in the basis set then these signals will not be 
periodic in the total sample interval. These signals which are not 
commensurate with the natural period exhibit discontinuities at their 
periodic extension. These discontinuities contribute to the spectral 
content of the signal calculated by the Fourier transform. The phenom- 
enon just described is a particular case of Gibbs phenomenon (Bracewell 
1978) , which refers to the sidelobe or ripple effect in the function 
produced by the discontinuity in the transform. It is this behavior in 
the transform, due to truncation of the autocorrelation, which prompts 
us to perform windowing. 

Because of the reasons stated above, Blackman and Tukey (1958) 
developed a technique whereby an autocorrelation function is modified by 
multiplying a weighting function in order to obtain a Power Spectral 
Estimate. This weighting function, which multiplies the autocorrelation 
in the function domain, is equivalent to the convolution of the Fourier 
transforms of the weighting function and the autocorrelation in the 
transform domain. The Fourier transform of the weighting function is 
called a window. The expression for the above relationship is the 
following: 



)exp(-i2^rsx)dx. 


( 16 ) 
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To avoid a convolution in the power spectral estimate calculation, 
Blackman and Tukey presented the method for windowing in the continuous 
case by multiplying the autocorrelation with the weighting function. 
The product of the above was then Fourier transformed to obtain the 
power spectral estimate. For the discrete case the convolution in the 
transform domain is not difficult. Therefore convolution of the window 
with the transform of the autocorrelation is an alternative. Robinson 
(1980) also elaborates on the above method. The specific methodology 
used in this work will be reviewed later. 

A number of windows have been introduced in the literature 
(Geckinli and Yavuz 1983). Harris (1978) elaborates on several of these 
and presents quantitative measurements which measure the performance of 
each. Table 1 of Harris's paper lists these parameters. As with the 
truncated input autocorrelation, the windows themselves will exhibit 
spectral leakage. This leakage is measured in Harris's work by the peak 
sidelobe level and the asymptotic rate of falloff of the sidelobes. 
These quantities are illustrated in Harris's Figure 5. Though the 
ratings of these quantities given by Harris are a measure of the per- 
formance quality of the window, they may not accurately reflect the 
usefulness of the windows in this work. The sidelobes play a critical 
role in the work presented here, because through these we will be 
determining the proper combinations to clear channels in the power 


spectral estimate. 


13 


CHAPTER 3 

CHANNEL CLEARING: THEORETICAL BACKGROUND AND INITIAL RESULTS 

In 1985 Smith introduced the technique of channel clearing. This 
method is based on the ability to apply weighting functions alone or in 
linear combinations to the autocorrelation function in order to clear 
channels or domain segments in the power spectral estimate. 

A specific problem as stated by Smith is the location of a small 
spike near a large spike in the power spectrum. In this case applica- 
tion of a window along with its sidelobes effect on the power spectrum 
has masked the smaller spike. Figure 3.1 shows the spectrum with a 
spike in channel 16 and a spike of l/20th amplitude in channel 20. 
Figure 3.2 shows the same spectrum with a single window applied. The 
influence of the larger spike combined with the window has masked the 
smaller spike in channel 20. The idea behind the technique is to 
combine windows linearly and apply to the spectrum containing the large 
spike only in order to clear or bring to zero the channel in which the 
smaller spike resides. When this combination is established, it is 
applied to the spectrum which contains both spikes. Subsequent to this 
the smaller spike would be detected. 

This combination of windows is called a hybrid window. We are 
allowed to make this construction because scaler addition and multipli- 
cation by a constant are preserved through the Fourier transform. This 
is shown in the following relationship 
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Z a . w i «i <»). 

where a is a coefficient such that the sum of all a's for a particular 
hybrid window equals unity. This normalization is recommended due to 
the fact that the zero lag of the autocorrelation is a maximum and 
amplification of this value can occur if the weighting is nor normal- 
ized. Another result that is related to this normalization is that the 
hybrid window w (x) is bounded by the windows of which it is composed. 
A consequence of this bounding is that it prevents amplification in any 
of the sidelobes . The proof of this bounding is given in Smith's thesis 
(1985). This proof extends to the case where the hybrid window is 
constructed of n component windows. In this case an envelope is con- 
structed with an upper limit window comprised of greatest components of 
all the input windows and a lower limit window comprised of the least 
magnitude components of the input windows. The analogy here results in 

the following statement: that for any x there are two component windows 

such that 


w,(x)< w fc (x)<* u (x) 

where u - the upper bound window, 1 - the lower bound window and 

h 

2>u>-i 

i'l 

and 


0< a(i)< 1. 
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The last property discussed in this section is the combination 
process and the use of hybrid windows as component windows in subsequent 
hybrid windows. It was shown that a linear combination of any number of 
windows, with normalized coefficients, produces a hybrid window. This 
hybrid window can be used in the construction of a new hybrid window, 
again with proper normalization of the coefficients. This property is 
important in the calculation of hybrid windows when attempting to clear 
more than one channel simultaneously. The technique used in determining 
the coefficients for the windows will demonstrate this. 


The initial results of this technique showed that a single channel 
c,n ba cleared within 10' 7 of the noraallzad peak value. The use of 
double precision may result in clearing within lo' 15 . The experiment 
was conducted by putting a unit spike in channel 31 and a spike of 
l/20th amplitude in channel 35. The purpose of the experiment was to 
see if a combination of windows could be calculated to eliminate the 
effects of the spike in channel 31 at channel 35. A combination of 
three windows were used in order to eliminate channel 35 from the 
spectrum with only the spike in channel 31 present. Windows 3. 6 and 9 
of Table 4.1 with respective coefficients of .0000016, .0000074, 

.9999910 were the choices. This hybrid window was then applied to the 
spectrum with the two spikes and clearly the spike in channel 35 was 
detected (Figure 5.3 and 5.4). The exact procedure for calculating the 
coefficients will be detailed in the next chapter. 

The second experiment involved replacing the spikes with Gaussians 
centered at channels 31 and 35. Again the maximum amplitude of each 
was 1.0 and .05 respectively. The purpose of this experiment is to 
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attempt to clear a finite extent or more than one channel and to treat 
peaks which are not spikes in the power spectral estimate. In this 
experiment a combination of two hybrid window were used to clear chan- 
nels 35 and 36. The final hybrid window was constructed from windows 2, 
6 and 9 from Table 4.1. The respective coefficients were .3533488, 
.0501631, and .5965287. The clearing of channels 35 and 36 was suc- 
cessful (Figure 5.6). A similar result is achieved in this work and is 
shown in chapter 5 . 
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CHAPTER 4 

METHOD AND ALGORITHMS FOR CALCULATIONS 

The method, at this time, is based on the assumption that we have 
some knowledge of the type of power spectral estimate and that portion 
relative to a large peak in which we need to consider the effects of 
windowing. For example, in the results reviewed in chapter 3 we set up 
the spectrum such that we knew channel 31 contained a dominant peak. We 
can model this response independently and design the appropriate hybrid 
windows to clear certain channels relative to the peak from the effect 
of this response. 

The exact method employed here utilizes a combination of hand and 
computer calculation. The first step in the method involves creating 
the model spectrum. For example, in the first experiments we placed a 
unit spike in a channel. Upon the creation of the model input each 
window is scaled to unity and applied to the model data. This is done 
by a program of the type listed in Appendix 1. These windows are listed 
in Table 4.1 and graphically represented in Figures 4.1 thru 4.11. The 
numeric output is listed for each window result. The actual calculation 
by the algorithm is as follows: 

1) Input window number (based on Table 4.1) 

2) Input of coefficients 

3) Algorithm multiplies windows by coefficients 

4) Algorithm inverse Fourier Transform of input 

5) Algorithm multiplies windows by inverse transformed 


input 
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6) Algorithm Fourier transform result to frequency domain. 

The output represents the numeric value for each frequency domain 
channel. From this a certain channel will be selected to extinguish 
This selection requires for the window combination that there is a 
negative and a positive result in that particular channel. This re- 
striction is due to the fact we are limiting the value of the coeffi- 
cients to being positive and less than unity (Smith 1985) . With the 
window and the associated values in the specific channel chosen we can 
use the following formula to calculate the coefficients. 


w(hybrid)-(l/(a(l)+a(2)))(a(l)W(2)+a(2)W(l)) 


(17) 


where 


(l/(a(l)+a(2)) ) -normalization. 

Therefore the coefficient applied to window 
W(l)-a(2)/(a(l)+a(2)) f and for 
W(2)-a(l)/(a(l)+a(2)). 

The following is a numeric example of this calculation. Suppose from 
the initial output we have values for channel 35 of -.0004536 and 
.0658214 from windows 1 and 2 respectively. Then in the formula 
a( 1)— . 0004536 and a(2)— .0658214. These values can be inserted into 
equation 17 . The normalization factor is equal to 
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l/( .0004536+. 0658214)-15. 0886457. 

Therefore the coefficient for window 1 is equal to 
a(2 ) 15 . 0886457- ( . 0004536) (15 . 0886457)- . 0068442 

and for window 2 is equal to 

a( 1) 15. 0886457-( .0658214) (15. 0886457)-. 9931558. 

These coefficients can then be applied to a spectrum with a response in 
the channel cleared other than the response resulting from the model 
input. This hybrid window along with another hybrid window clearing the 
same channel may be utilized in the construction of new hybrid windows 
to extinguish further channels. We are able to do this because we know 
that for the first channel cleared there are zeroes or relatively small 
numbers resident in this channel. So when we add and then apply by 
multiplying these hybrid windows, the first channel will remain extin- 
guished. The same process for clearing as outlined above, in which we 
can utilize the numeric output from the first hybrids to identify a 
channel with a positive and negative value, is used. 
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CHAPTER 5 


COMPUTATIONAL RESULTS 

The first experiment conducted was to duplicate the results 
obtained by Smith (1985) in his work, which was reviewed in Chapter 3. 
The purpose for this exercise was twofold. 1) The algorithm utilized 
by Smith was transported from a mainframe environment to a personal 
computer. The Fortran used by the two systems is different. As a 
consequence of this, some conversion was necessary. 2) The plots 
produced in the original work displayed the resultant Power Spectral 
Estimate with the normalized amplitude values from the calculation. In 
keeping with the literature, a dB-down representation of the results is 
generated for this work. We wanted to present Smith's results in this 
manner also. 

Figure 5-1 duplicates the result of the window combination used to 
extinguish channel 35 with the unit spike in channel 31 only. Notice 
from this scale the extinction of channel 35 is not observable. The 
rescaling to dB-down (Figure 5.2) clearly shows the extinction at 
channel 35. Figures 5.3 and 5.4 show the results of the applied hybrid 
window to the spectrum consisting of spikes in both channels 31 and 35. 
The spike located in channel 35 is l/20th the amplitude of the spike in 
channel 31. Restating the conclusions from the original results, it 
appears that the spike in channel 35 suffered no adverse effects gener- 
ated by the spike in channel 31 subsequent to windowing. 
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Figures 5.5 through 5.8 duplicate the results for the Gaussian 
input introduced in chapter 3. The numeric objective in this experiment 
was to clear channels 35 and 36. Again in the amplitude plot of Figure 
5.5, the extinction of these channels is not detectable. In Figure 5.6 
channels 35 and 36 are obviously down from their adjacent channels. 
Figures 5.7 and 5.8 show the results of application of the same hybrid 
window to the two Gaussian input. Notice the increase of amplitude in 
channel 34 when comparing Figures 5.6 and 5.8. This increase in ampli- 
tude is attributed to the application of the hybrid window to clear 35 
and 36 but not 34 and presence of the smaller Gaussian centered in 
channel 35. 

As explained in the methodology in chapter 4, it is necessary to 
examine the actual numeric output in order to calculate the coefficients 
for the windows. In the subsequent experiments these lists of values 
will be given as tables. 

The main purpose of the experimental work is to introduce new data 
input to this technique in order to examine the robustness thereof. 

first computation is to investigate, for the same magnitude 
input spikes used in the above, whether the channel clearing is affected 
by the location of the spikes in the spectrum. The first data set was 
constructed by placing the large spike in channel 16 and the smaller 
spike in channel 20. Application of the same scaled windows were 
applied to these data with the intent of clearing channel 20 from the 
effects of the spike in channel 16. The results are displayed in 
figures 5.9 through 5.12. From analyzing figure 5.10 it is evident that 
the extinction at channel 20 is not significant relative to its two 
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adjacent channels. This is borne out in the output run listed in table 
5.1. Although the extinction was not as complete as it was in channel 
31-35 case, the sidelobe level of the spike in channel 16 generated by 
the windowing did not mask the smaller spike in channel 20 (Figure 5.11 
and 5.12). The technique was then applied to these same input data with 
windows 3 and 9 applied with respective coefficients of .000002279 and 
.999997736. Table 5.2 lists the results of the application of these 
windows to the input of a spike in channel 16 only. Comparing channels 

19 through 21 in tables 5.1 and 5.2 demonstrates that the new calculated 
windows performed better. Figures 5.13 and 5.14 also show that channel 

20 has been better cleared relative to its adjacent channels, but once 
again, it does not appear that the sidelobe effect of the spike in 
channel 16 is greatly influential out to channel 20. Because of this 
the application of the windows to the input of spikes in both channels 
16 and 20 was not performed, knowing the smaller spike in 20 would be 
detectable . 

The second experiment involves taking different percentages of 
possible lag values of the autocorrelation as represented in the input 
power spectral estimate. By placing spikes in the spectrum, as we have 
in the first experiments, we have represented an autocorrelation func- 
tion from a very small percentage of total possible lags (re. Chapter 
2). Typically in the use and application of the autocorrelation func- 
tion a greater percentage of lags are used to represent the function. 

In this work we will examine the effect of taking 20 and 100 percent of 
the lags. The effect of taking different percentages of lags in itself 
is a type of windowing. Figures 5.15 and 5.16 show the function domain 
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weightings representing 20 and 100 percent lags respectively. These 
weightings are multiplied times the input similar to the other function 
domain windows to give the effect of using 20% or 100% of available 
lags . 

The first analysis was performed with a single spike in channel 31 
using 20 percent of the lags. Again the windows will be designed to 
extinguish channel 35. Figures 5.17 and 5.18 show the smoothing of the 
spectrum from a larger number lags. It should be noted here that, when 
adding the smaller spike to the spectrum, the smoothing effect of taking 
20 percent of the lags is present as a result of the presence of the 
smaller spike also. Consequently, there is a larger effect in channel 
35 (Table 5.3). The clearing of channel 35 was accomplished twice using 
two different combinations of windows. At this point, for the remaining 
analysis, the dB-down representations will only be shown with the 
exception of Figure 5.24. The first result used windows 1 and 6 (Table 
4.1) with respective coefficients of .312745355 and .687254645. Figure 
5.19 and Table 5.14 give the results of this calculation. The extinc- 
tion of channel 35 was accomplished along with channel 27. It should be 
stated that the calculations were not intended to clear channel 27, but 
the clearing appears to be the result of symmetry. Figure 5.20 shows 
the result of applying the above window combination with the input 
consisting of the smaller spike in channel 35. The spike in channel 35 
is now detectable apparently with the correct relative amplitude. The 
second combination of windows applied to this data set was 3 and 7 with 
respective coefficients of .450843611 and .549156389. Figures 5.21 and 
5.22 along with Table 5.5 show a similar result as for window 
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combination 1 and 6 . The last part of this experiment was to attempt to 
clear a channel closer to the main peak. The input is the 20 percent 
lags of the autocorrelation with the spectrum containing the large spike 
in channel 31 and the smaller spike in channel 35. In the previous part 
of this experiment two window combinations were given to clear channel 
35. In order to clear the second channel the method of chapter 4 will 
be used in conjunction with the fact that the combination of existing 
hybrid windows can be used as a new hybrid window. An important aspect 
of methodology stems from this analysis and will be shown in detail in 
another experiment. It is important when calculating the coefficients 
for the extinction of the second channel that the presence of the other 
spike, which may be detected due to the clearing of the first channel, 
be taken into account. The two hybrid windows used consisted of windows 
1, 6, 3 and 7 with coefficients .1325286, .291230213, .259794658, and 

.316446529 respectively. Figure 5.23 and Table 5.6 show that this 
combination cleared channel 33 and channel 35, resulting in the preser- 
vation of the small spike. 

The second experiment involves the use of 100 percent of the lags 
of the input data. This input is shown in Figures 5.24, 5.25, and Table 
5.7. The smoothing effect from using all the lags is greater. Because 
this broadening of the main lobe, this data set was used more exten- 
sively in this analysis. It is felt that this broadening effect can 
better demonstrate the robustness of the method. The experiments 
performed this on this data set required the use of many different 
window combinations. Therefore, in like manner to the 20 percent lag 
case, several window combinations were used to clear a single channel. 
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Table 5.8 gives a list of window combinations and their respective 
coefficients which clear channel 35. This table also provides the 
Figure and Table numbers associated with the results of these calcula- 
tions. The five window combinations utilized were all successful in 
clearing channel 35 to at least -70 db. The first four combinations 
show symmetry in their spectra. The fifth window combination is anoma- 
lous in this regard. This variability in symmetry should be investi- 
gated but was not in this analysis. 

The second a part of this analysis was to clear two consecutive 
channels or in effect clear a small domain segment. For this we uti- 
lized the window combinations of (3,8) and (7,10). The coefficients are 
as follows: a(3)-. 201309972 , a(8)-. 193117857 , a(7)-. 494649063 , 
a(10). 110862395. Figure 5.31 and Table 5.14 demonstrate the successful 
results of this combination. The concept of clearing domain segments or 
successive channels has more importance in considering broad lobed input 
spectra (e.g. broad Gaussians, Smith 1985). 

The third part of this analysis refers back to the first experi- 
ment. It demonstrates the importance of the methodology of the tech- 
nique when clearing more than one channel after detecting an event in 
the first channel cleared. The combination of hybrid windows (1,6) and 
(4,5) will be used to clear channel 44 after clearing of channel 35. 
Figure 5.32 and Table 5.15 show the extinction of channels 35 and 44 in 
the presence of a spike in channel 31 with 100 percent of lags. The 
following coefficients were used: a(l)- . 517188658 , a(6)- . 007515903 , 
a(4)-. 307604029, a(5)-. 16769141 . Figure 5.33 and Table 5.16 show the 
results of using the same windows and coefficients applied to the input 
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spectrum including the smaller spike in channel 35. The presence of 
this spike has adversely affected the ability of these specific coeffi- 
cients to clear channel 44. Recalculation of the coefficients, utiliz- 
ing the same window combinations, produces the results shown in Figure 
5.34 and Table 5.17. These recalculated coefficients are a(l)-. 5290774 , 
a(6)-. 007688673, a(4)-. 29979798 , a(5)-. 163435915 . This procedure has 

more importance as the second channel to be cleared comes nearer to the 
secondary peak. Also, the magnitude of the secondary peak influences 
this result proportionally. The next experiment extends this technique 
to clearing three channels simultaneously without the presence of a 
secondary or tertiary response in any of the previously cleared chan- 
nels. The hybrid window combinations of (1,6) (4,5) (7,10) are 

used in this case with the following coefficients: 

a(l)-. 661751704 
a(6)-. 009616726 
a(4)-. 198866131 
a(5)-. 108412565 
a(7)-. 017443401 
a(10)-. 003909473 

Figure 5.35 and Table 5.18 show the success of clearing channels 35, 
39, and 44. Any further extension of this to more channels or larger 
domain segments will necessitate the automation of the method. This is 
due to the number of computations to be made. 
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The last experiment considers three large spikes in the input 
spectrum at channels 30, 35, and 40. In this case we considered a small 
percentage of lags. Using windows 3 and 8 with respective coefficients 
.785292479 and .214707521, the attempt to clear channel 33 was made. 
Figure 5.36 and Table 5.19 show the extinction of channel 33 greater 
than -70 dB. Figure 5.37 and Table 5.20 show the results of the same 
window combination with a l/20th amplitude spike placed in channel 33. 
The smaller spike is detected despite the proximity of the two large 
spikes . 
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CHAPTER 6 

CONCLUSIONS AND FUTURE CONSIDERATIONS 

The purpose for and technique of channel clearing, in the power 
spectral estimate by application of linearly combined window functions, 
has been reviewed. It has been established in previous work, that this 
method is effective in clearing up to two channels. This work extends 
the applicability of this method to different input data. 

The first calculation demonstrated that the location of the spike 
in the input spectrum varied the result of channel clearing calculations 
relative to an equivalent valued spike used in the previous work. Part 
of this effect could be the result of the transform wrap-around. 

The second investigation involved the use of varied percentages of 
lags relative to the input autocorrelation and its associated power 
spectral input. It was demonstrated that the effect of broadening the 
main lobe did not effect the channel clearing capabilities of the 
method. Several calculations were performed using the greater percent- 
age lag input. It was shown that three channels could be cleared 
simultaneously. Two consecutive channels or a small domain segment were 
cleared. The importance of the methodology in the calculations was 
demonstrated by clearing two channels simultaneously with and without 
the presence of a smaller spike in the first channel cleared. The final 
calculation demonstrated the effectiveness of the technique when three 
large spikes are relatively close together. A single channel between 
two of the large spikes was cleared. 
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The method of channel clearing in the power spectral estimate needs 
further evaluation. Specific to this is the need to examine varying in- 
put data, different windows, more arithmetic precision and further 
clearing of more channels or larger domain segments. In light of the 
latter it is important that the technique become more automated. This 
automation could include the capability of channel marching, whereby 
each channel is cleared successively. Lastly, the implementation of 
this technique to solve real data situations needs to be investigated. 
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TABLE 4.1 


WINDOW NUMBER 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


w(x) for -1/2 <x<+1/2 
1 

l-4x 2 
cos (irx) 

Sin(27rx)/27rx 

1+2 | x | 

.5+.5Cos(2wx) 

. 54+ . 46Cos ( 2*77 x) 

1-2|x|Cos(29tx) +(l/7r) | Sin(2?rx) | 
.42+. 5Cos(2wx)+. 08Cos(4ttx) 

1 ~ 2 1 x | 

. 5+. 5Cos(2w( | x | - . 15) )/ . 7 
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TABLE 5. 1 


CHANNEL 

AMPLITUDE 

1 

0 • 0000000 1 6B 

Am 

0. 0000000108 

T, 

0.0000000044 

4 

-0. 0000000199 

5 

0 . 0000000009 

6 

0. 0000000029 

7 

-0 . 0000000 1 7 

e 

— 0. 0000000442 

9 

— 0. 0000000046 

10 

-0 . 0000000 176 

1 1 

0. 0000000126 

12 

-0. 0000000458 

13 

0. 0000000537 

14 

0. 079999164 

15 

0. 49999979 

16 

0. 8400017 

17 

0. 49999982 

18 

0.079999164 

19 

0 . 000000 1 0 1 2 

20 

0 • 0000000 1 84 

21 

0. 0000000308 

n 

A. Am 

— 0 . 0000000 167 


0 • 0000000334 

24 

0 . 000000053 1 

25 

0 . 000000020 1 

26 

0. 0000000066 

27 

0. 0000000143 

28 

0 . 0000000 1 07 

29 

0. 0000000075 

30 

-0. 0000000085 

31 

-0. 0000000024 

32 

-0 . 0000000034 


CHANNEL 

AMPLITUDE 


0 . 0000000 1 3 

34 

0. 0000000086 

35 

0 . 0000000253 

36 

0. 0000000008 

37 

0 . 000000006 1 

38 

0 . 0000000 1 38 

39 

0 . OOOOOOO 143 

40 

0.0000000207 

41 

-0 . 000000007 4 

42 

0.0000000047 

43 

0 . 0000000039 

44 

0 • 0000000 121 

45 

0. UOOOOU0O4B 

46 

0. 0000000053 

47 

0 . 0000000 1 33 

48 

—0. 0000000037 

49 

0.0000000298 

50 

0. 0000000596 

51 

O. 0000000447 

52 

0. 0000000075 

53 

— O. 0000000034 

54 

0. 0000000006 

55 

0 . OOOOOOO 173 

56 

0. 0000000086 

57 

0 . 0000000068 

58 

0 . OOOOOOO 1 36 

59 

0 . OOOOOOO 1 02 

60 

0. 0000000103 

61 

0 . OOOOOOO 1 67 

62 

-0 . OOOOOOO 1 02 

63 

0.0000000017 

64 

-0 . 000000020 1 
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TABLE 


■NNEL 

AMPLITUDE 

CHANNEL- 

1 

0. 00000uuO4 2 



—0. 0oooO<. # 0O-> 

34 

3 

0 • 0000 1 -*00 1 

35 

4 

—0 . 0000000 158 

36 

5 

9 . 3563823E- 1 1 

37 

6 

—0 . 0000000035 

38 

/ 

—0. 0000000 1 24 

39 

8 

-0. 000 0000 ©4 3 

40 

9 

0 • 0000000' '6 

41 

10 

— 0 . OOOOOuO 174 

42 

1 1 

0 • 0000000309 

43 

12 

—0 • 0000'.*006 

44 

13 

0. 000O000731 

4 5 

14 

0. 079999643 

46 

15 

U. 49999976 

47 

16 

0 . 6400008 1 

48 

17 

O. 4999997o 

49 

18 

0.079999611 

50 

19 

0 • 000000 1 086 

51 

20 

0 . 0 0 0 0 0 0 00 4 2 

52 

21 

0. 0 00 0 00 0 399 

53 

22 

— o . 0000000 152 

54 


0 . 00000004 1 5 

55 

24 

O. 0000000408 

56 

23 

0 . OOOOOOO 2 U 4 

57 

26 

0. 0000000044 

58 

27 

0 • OOOOOOO 1 6 

59 

28 

0. 000000O065 

60 

29 

-0. 0000000057 

61 

30 

-0 . OOOOOOO 1 1 

62 

31 

-0.0000000135 

63 

32 

0. 0000000003 

64 


AMPLITUDE 

0 . 0000000033 
O. 00000001 19 
0. 0000000155 
■0 . 0000000069 
0 . 0000000 1 32 
0 . 00000001 19 
0.00000002 1 6 
0 . 0000000 1 32 
0. 0000000066 
-0. 000000004 
0. 0000000079 
0 . 00000001 19 
0. 0000000102 
0. 0000000046 
0. 0000^00123 
0 

0. 0000000 29B 
0 

0 . 0 00 Oo 00 299 
0. 0000000 1 96 
O. 000000002 1 
0. 0000000069 
0. 0000000212 
0. 0000000042 
O • 0000000 117 
0. 0000000097 
— 0 • OOOOOOOO 1 3 
O . 0000000 147 
O. 0000000143 
O . OOOOOOOO 1 9 
-0.0000000091 
-0.0000000174 
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TABLE 


CHANNEL 

AMPLITUDE 

1 

0 

2 

o 

7 . 

0 

4 

0 

5 

0 

6 

o 

“7 

/ 

o 

8 

o 

9 

0 

10 

0 

1 1 

0 

12 

0 

13 

0 

14 

0 

15 

0 

16 

o 

17 

0 

IB 

0 

19 

0.0032585 

20 

0.036195 

21 

0. 0032575 

no 

< 4 . 4 . 

0 • 0037895 

23 

0 . 0032575 

24 

0. 0041285 

25 

0- 00325S5 

26 

0.0049645 

27 

0.0032575 

28 

0 • 0080555 

29 

0. 0032585 

30 

0. 0473235 

31 

1 

32 

0. 0473235 


CHANNEL 

AMPLITUDE 

33 

0.0032585 

34 

0 . 0080555 

35 

0. 0032575 

36 

0.0049645 

37 

0. 0032585 

38 

0.0041285 

39 

0.0032575 

40 

0. 0037895 

41 

0. 0032575 

42 

0. 0036195 

43 

0. 0032585 

44 

0 

45 

0 

46 

0 

47 

0 

48 

0 

49 

0 

50 

0 

51 

0 

52 

0 

• 53 

0 

54 

0 

55 

o 

56 

0 

57 

o 

58 

0 

59 

o 

60 

o 

61 

0 

62 

0 

63 

o 

64 

0 


[>iG£ IS 

OF POOR QUALITY 
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.ABLE 5.4 



CHANNEL 

AMPLITUDE 

CHANNEL 

AMPLITUDE 

— 

1 

-0.0087503614 

33 

0.014556956 


xl 

0 . 0087503623 

34 

0.02156426 


“7 

—0 . 00875032 1 3 

35 

0. 0000000596 


4 

0. 0087503623 

36 

0.017506585 


5 

-0.0087503539 

37 

-0.0013481798 


6 

0. 008750366 

38 

0.016409133 


7 

— 0. 0087503577 

39 

-0.0017531442 

— 

8 

0. 0087503577 

40 

0.015963767 


9 

-0. 0087503418 

41 

-0. 0019281537 


10 

0. 0087503688 

42 

0.015740924 

_ 

1 1 

—0. 0087503558 

43 

-0. 0032289964 


12 

O. 0087503623 

44 

0. 0098700784 


13 

—0. 00875037 16 

45 

-0. 008750351 1 


14 

0. 0087503595 

46 

0.0087503549 


15 

-0. 0087503456 

47 

-0.0087503772 


16 

0.0087503605 

48 

0.0087503623 


17 

—0. 0087503577 

49 

— O. 0087503707 

— 

18 

0. 0098700728 

50 

0. 0067503605 


19 

-0.0032290181 

51 

-0. 0087503623 


20 

0. 015740927 

52 

0.0087503623 


21 

-0. 001928166 

53 

-0 . 00875036 1 4 


-v*> 

X X. 

0. 015963722 

54 

O . 008750352 1 


2 •-» 

—0.001 7533107 

55 

— O. 006750338 1 


24 

0.01 6409073 

56 

0. 0087503679 

— 

nr. 

-0.0013481884 

57 

-0. 0067503623 


26 

0. 017506564 

58 

0.0087503521 


27 

—0. 0000000758 

59 

-0.0087503642 


28 

0.021564264 

60 

0.0087503688 


29 

0. 014556956 

61 

-0.0087503502 


30 

0.41562104 

62 

0. 0087503726 


31 

1.336518 

63 

-0. 0087503688 


32 

0.41562104 

64 

O. 0087503688 


C'f 


QUALITY 
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TABLE 5.5 


CHANNEL 

AMPLITUDE 

CHANNEL 

1 

-0. 0015726413 

33 

2 

0. 00 1 5736328 

34 

T ( 

—0. 0015765488 

35 

4 

0. 00 1 58 1 52 1 5 

36 

5 

—0 . 00 1 58B5B 1 . >3 

37 

6 

0 . 00 1 59S0 1 69 

38 

"7 

-0. OOl 6i.)99684 

39 

B 

0. 0016247729 

40 

9 

-0. 0016428066 

41 

10 

0 . 00 1 6647 036 

42 

1 1 

-0. 001691038'/ 

43 

12 

0. 0017228881 

44 

1 3 

-0. 0017615802 

45 

14 

0. 00 1 609 1 003 

46 

15 

—0 . 00 1 869089 

47 

16 

<-> . oO 19 4 7 6o7 / 

4S 

17 

-0. 0020815721 

49 

18 

0.0034399191 

50 

19 

0 . 003 1 004637 

51 

20 

0. 0093984222 

5^ 

21 

0.0042746561 

53 


0. 010015843 

54 

— N T 

4L ■-> 

0. 0038203243 

55 

24 

0.01 1394353 

56 

25 

0.0028025657 

57 

26 

0 . 0 1 48 1 205 

58 

27 

-0. 0000031461 

59 

28 

0.027611868 

60 

29 

-0.010478545 

61 

30 

0.49984768 

62 

31 

1 . 2078092 

63 

"T 

0.49984694 

64 


AMPLITUDE 
-O. 010477016 
O. 027609587 
O 

O. O 1 4808223 
0.0026072514 
0.011388922 
0.0038268473 
0. 010008626 
O. 004282867 1 
0. 0093892235 
0.0031106449 
0.0034286149 
-0.0020691007 
0.0019359596 
-0.001854112 
0. 0017926926 
— O . 00 1 7 4 3 c 6 2 2 
O. 0017033712 
-0.0016697475 
0.0016414967 
-0.0016175528 
0. 001597214 
-0.0015799081 
0.0015652351 
-0.0015527973 
0.0015423266 
-O. 001533656 
0.0015265713 
-O. 00152091 
0. 0015166529 
-0.0015136404 
0.0015118601 
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TABLE 5.6 


CHAN! 


EL 

AH PL I TUDE 

CHANNEL 

AMPLITUDE 

1 

-0. 0046152702 

33 

-0 . 00000002 1 


0.0048156313 

34 

O . «.>456839o5 

-r 

—0 . 0048 1 752o7 

3 b 

0. 055507839 

4 

0. 0048205 

36 

0.036541209 

5 

— 0. 00482467 89 

37 

0. 001 1355346 

6 

0 • 0048302528 

36 

0. 014625136 


-0. 0048373439 

39 

0.0014689919 

8 

0. 004846076 

40 

0.013229766 

9 

-U. 0048567597 

41 

0.001fc927826 

10 

0.0048696799 

42 

0. 01266754 

1 1 

— 0 . 0048852339 

43 

0.0004212931 

12 

0. 0049040862 

44 

0. 00644B7387 

13 

-0. 0049268967 

4b 

-0.0051310528 

1 4 

0.0049549919 

46 

0.0050476063 

15 

-0. 0049903756 

47 

-0.0049959356 

16 

0 . 0050380644 

48 

0.0049571553 

17 

-0.0051 163761 

49 

-0.0049262699 

18 

0 . 006427 1176 

50 

0. 0049009654 

19 

0. 0004415293 

51 

-0. 0048799003 

20 

0.012593273 

52 

0.0048622494 

21 

0.0017247929 

53 

-0.0048473131 

'jO 

0.013049007 

54 

0.0048346845 


0. 00 1 534 2634 

55 

—0. 0048238523 

24 

0 . 0 1 4043696 

56 

0.0048148078 

2 b 

0. 001 1 129705 

57 

—0 . 004807 068 

26 

0.016499601 

56 

0.0048006219 

27 

0. 0000582 112 

59 

-0.0047951946 

2B 

0.025635593 

60 

0.0047908537 

29 

0 . 000 171 23 1 2 

61 

-0.0047873775 

30 

0. 46460469 

62 

0.0047847191 

51 

1 . 2619438 

63 

-0.0047829049 


0. 46447643 

64 

0.0047817547 


TABLE 5.7 


CHANNEL AMPLITUDE CHANNEL AMPLITUDE 


1 

0 


o 


0 

4 

0 

5 

0 

6 

0 

7 

0 

e 

0 

9 

0 

10 

0 

1 1 

0 

12 

o 

13 

0 

14 

0 

15 

0 

16 

<7 

17 

0 

18 

0 

19 

0. 0000972 

20 

0 • 0033475 

21 

0. 0000965 


0 . 0050 195 

' v T 

X--> 

0. 0000955 

24 

0. 0083305 

25 

0. 0000945 

26 

0. 0163905 

27 

0 . 0000935 

28 

0 - 0456595 

29 

0 „ 00009 1 5 

30 

0.4115925 

31 

1 


0.4115925 



0 . 00009 1 5 

34 

0.0456595 

35 

0. 0000935 

36 

0. 0163905 

37 

0. 0000945 

38 

0. 0083305 

39 

0. 0000955 

40 

0. 0050195 

41 

0. 0000965 

42 

0. 0033475 

43 

0. 0000975 

44 

0 

45 

0 

46 

0 

47 

0 

48 

0 

49 

o 

50 

0 

51 

0 

52 

0 

53 

0 

54 

o 

55 

0 

56 

0 

57 

0 

58 

o 

59 

0 

60 

0 

61 

0 

62 

0 

63 

0 

64 

o 
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LU 
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TABLE 5.9 



CHANNEL 

AMPLITUDE 

CHANNEL 


1 

- 0.000630172 

33 



0.000630194 

34 



- 0. 0006300242 

35 

— 

4 

0 . 000630 1 73 1 

36 


5 

- 0 . 000630 1 632 

37 


6 

0 . 0006301354 

38 


7 

- 0. 0006301847 

39 


8 

0. 0006301642 

40 


9 

- 0 . 000630 1 1 93 

41 


10 

0 . 0006301703 

42 


1 1 

- 0. 0006301042 

43 


12 

0. 0006301827 

44 


13 

- 0 . 00063-01479 

45 

— 

14 

0. 0006301561 

46 


15 

- 0 . 000630 1 1 06 

47 


16 

0. 0006301499 

48 


17 

- 0. 000630165 

49 


18 

0 . 0006308669 

50 


19 

- 0 . 0004 1 25 1 3 

51 


20 

0.0072785877 

52 

— 

21 

— 0.0003766093 

53 


2 . *L 

0.010598592 

54 


23 

— 0. 0 0 03450355 

55 

L 

24 

0 .0171 73 1 54 

56 


25 

- 0. 0002654459 

57 


26 

0.03317773 

58 


27 

- 0. 0000001369 

59 


28 

0 . 09 1 296434 

60 


29 

0.0028264024 

61 


30 

0. 82508206 

62 

— 

31 

1 . 99094 1 

63 


32 

0. 82508206 

64 


- ' •& 




Y 


AMPLITUDE 
0.0028264243 
0.091296464 
•0.0000001788 
0.033177733 
0.0002654514 
0.017173227 
0. 0003448348 
0. 010598686 
■0. 0003786089 
0. 00727861 
■0.0004125431 
O. 0006308775 
■0. 0006301322 
0 • 000630 1564 
-0. 0006301673 
O . 000630 1616 
-0.0006301492 
O . 000630 1531 
-0 . 000630 1 033 
O . 000630 1 809 
-0.0006301597 
0 . 000630 1677 
-0 . 000630 1397 
O ■ 000630 1698 
-O . 000630 1646 
O. 0006301806 
-O . 0006 301446 
0.0006301526 
— O . 000630 1 35 
O . 000630 141 
-0.000630094 
O. 0006301437 


l~i in 


9 


AO 


FABLE 5.10 


CHANNEL 

AMPLITUDE 

CHANNEL 

1 

—0.000630 1 599 

33 


0 - 000630 1 826 

34 


— 0. 000630054 

35 

4 

0. 000630 1713 

36 

5 

-0. 0006301547 

37 

6 

0 « 000630 1 39 1 

36 

/ 

— 0 . 000630 1 939 

39 


AMPLITUDE 
0.0028263982 
0.09129646A 
— 0. OOOOOOl 192 
0.033177733 
-O. 0002654495 
0.017173231 
— O. 000344B291 


8 

0. 0006301592 

40 

0.01 0598685 


-0. 0006301 3 1 e 

41 

- 0 . 0 0 0 378608 2 

10 

0 . 000630 1 663 

42 

0. 0072786137 

1 1 

-0. 0006300929 

43 

-0. 0004125502 

12 

0. 0006301932 

44 

0. 0006308747 

13 

-0. 0006301431 

45 

-0 . 000630 1298 

14 

0. 0006301478 

46 

0 . 000630 1 54 

15 

-0. 0006301464 

47 

-0 . 000630 1718 

16 

0 . 000630 1 695 

46 

0. 0006301422 

17 

-0. 000630 149 

49 

—0 . 000630 1 503 

18 

0. 0006308607 

50 

0.0006301635 

19 

-0. 0004 1 25077 

51 

-0. 000630097 1 

20 

0. 00727858o3 

52 

0 . 000630 1 824 

21 

-0. 0003766138 

53 

-0 . 000630 1414 

O 

0.010598598 

54 

0. 0006301624 

“T 

-0. 0003450263 

55 

-0 . 0006301602 

24 

0.017173145 

56 

0 . 000630 1712 

^‘5 

-0. 0002654413 

57 

-0 . 000630 1574 

26 

0.033177722 

58 

0. 000630167 1 

27 

- 0. 0000003 42 1 

59 

-0.0006301 196 

20 

0. 09129e>44 9 

60 

O. 0006301573 

29 

0.0028263789 

61 

-0 . 000630 1 436 

30 

0. 82508206 

62 

0 . 000630 1451 

31 

1 . 99094 1 

63 

-0 . 000630 1 208 

32 

0 . 62508206 

64 

0 . 000630 1316 
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TABLE 5.11 


CHANNEL 

AMPLITUDE 

1 

-0. 00063341 96 

2 

0. 0000694477 

T 

-0.0008376982 

4 

0.0000603489 

5 

-0. 00085121 1 

6 

0. 000041 149 

7 

-0.0008749869 

8 

0. 0000099249 

9 

-0.00091 14152 

10 

—0 . 000037 1 338 

1 1 

-0. 0009647381 

12 

-0. 0001066394 

13 

-0. 0010416404 

14 

-0. 00021 10194 

15 

-0. 001 1546231 

16 

—0. 0003783755 

17 

-0.0013262827 

16 

-0. 000709017 1 

19 

—0. 0009096047 


0.0049861227 
— O . O C> O 6 0 O 9785 
0 . 007 1764561 
-0. 000664264 
0.01 167044 
-O. 0006962649 
0.022697762 
— 0. 0000034742 
0.0634931 18 
O. 030992692 
0. 88186336 
1 . 9435923 
0. 88186419 


CHANNEL 

AMPLITUDE 

33 

0. 030994404 

34 

0. 063495696 

35 

0 

36 

0. 022702038 

37 

-0 . 00069306 1 4 

38 

0.011676647 

39 

-0. 0006569907 

40 

0. 0071865907 

41 

-0. 0005918291 

42 

0. 0049963323 

43 

-0.0008982156 

44 

-0.0006964591 

45 

—0 . 00 131 24065 

46 

-0. 0003632295 

47 

-0.0011379927 

48 

-0. 0001928792 

49 

-0. 0010219279 

50 

-0. 0000849636 

51 

-0. 00094 1 0232 

52 

-0. 00001 1332 

53 

—0. 0008832853 

54 

0 . 0000406222 

55 

-0.0008414645 

56 

0.0000777696 

57 

-0. 00081 1311 

58 

0. 0001041 191 

59 

-0. 0007900024 

60 

0. 0001222219 

61 

-0. 0007757854 

62 

0 . 000 1 3358 1 1 

63 

-0. 0007676826 

64 

0. 0001390325 
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TABLE 5. 12 


EL 

AMPLITUDE 

CHANNEL 

i 

-0. 0006275484 

T "! 

***> 

0.0006295495 

34 


-0.0006271469 

35 

4 

0. 0006292728 

36 

5 

— 0 • 0006263803 

37 

6 

0. 0006283546 

38 

7 

-0. 0006246361 

39 

8 

0 . 000626748 

40 

9 

-0. 0006214 118 

41 

10 

0 . 000623729 1 

42 

1 1 

—0. 0006154042 

43 

12 

0. 000617 9126 

44 

13 

— 0. 0006037591 

45 

14 

0- 00060570 1 9 

46 

1 5 

— 0 . 0 y J 0 5 7 7 8 c- 1 J 7 

47 

16 

0. 0005774596 

48 

17 

-0. 000500306 1 

49 

18 

0. 000 56 6 5B 33 

50 

19 

—0 • 00034491 89 

51 

20 

0. 00733997 1 1 

52 

21 

— 0. 0003990874 

53 

27 

0. 0 1 075 148 

54 

2 

-0. 0003O98U6 

55 

24 

0. 017589815 

56 

■— > r- 
X w> 

-0.0001949752 

57 

26 

0.034976389 

58 

27 

—0 . 000000 1 007 

59 

28 

0. 10816382 

60 

29 

—0. 004 1 795578 

61 

30 

0. 80595249 

62 

31 

2 . 004 1 552 

63 


0.80595243 

64 


Al'IPL 1 TUBE 
— 0. 0041795401 
0. 106 1 6383 
O. 000000 1 1 92 
0. 03497e363 
-0. 0001949549 
0.017589815 
-0.0003095506 
O. 01075148 1 
- 0 . 0 00 399 0 322 
0.0073398859 
-0.0003448767 
0.0005664784 
-0.0005002396 
O. 00057731 17 
-O. 000577857 
0. 000605564 1 
—0. 0006037 396 
0 . 0006 1777 
-0.0006154663 
O • 0006 2 36 1 24 
-0 . 00062 1 4692' 
O. 0006266517 
— 0. 000624793 
0.00062831 1 8 
-0.0006267225 
0.0006292514 
-O. 0006278249 
0.0006297429 
-0.0006264985 
O. 000630028 
-0. 0006287869 
0. 000630223 1 
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TABLE 5. 13 


CHANNEL 

AMPLITUDE 

1 

- 0. 01612797 

r> 

— 0. 013706565 

7 

- 0.0089507475 

4 

- 0. 006047863 

5 

- 0.0060551067 

6 

- 0.012974271 

7 

- 0.016506637 

B 

- 0. 01492344 

9 

- 0. 0096571175 

10 

— 0 . 0052 1 1 600 1 

1 1 

- 0.0061681421 

12 

- 0 . 01 1894779 

13 

— 0 . 0 1 7562045 

14 

- 0.017429573 

15 

- 0.01 1094769 

16 

- 0.003568052 

17 

- 0.0021740822 

18 

- 0.0069588314 

19 

- 0. 017194647 

20 

- 0. 018569034 

21 

- 0.010338039 

X- 

0.0057766885 

•— > “5 

0. 015224451 

24 

0 . 01070234 1 

x~ J 

- 0. 0074134166 

26 

- 0.017665148 

27 

0. 0032037746 

26 

0.072982423 

29 

0. 29667765 

30 

1 . 0420661 

31 

1 . 616786 

32 

1 . 0462196 


CHANNEL 

AMPLITUDE 

33 

0. 30148342 

34 

0.07438153 

3S 

0 

36 

-0. 022786975 

37 

-0.01013279 

38 

O. 01270414 

39 

0. 020290535 

40 

0.0096400874 

41 

-0. 01097194 

42 

-0. 023223832 

43 

-0.021978475 

44 

-0. 0098132491 

45 

O. 001686646 

46 

0. 0017857002 

47 

—0 . 0087 6 1 9368 

48 

-0.020197596 

49 

-0.023196142 

50 

-0.015668012 

51 

-0. 0048406716 

52 

0.0002371874 

53 

-0.0046547656 

54 

-0. 014672683 

55 

-0. 021393139 

56 

-0.019007541 

57 

-0.010122581 

58 

-0. 0022366194 

59 

-0. 002312514 

60 

-0. 0098404614 

61 

-0.016509429 

62 

-0.020614322 

63 

-0.014695399 

64 

-0.0056453291 
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1 ABLE 5.14 


CHANNEL 

AI1PL I TUDE 

CHANNEL 

AMPLITUDE 

1 

-0.010013235 

33 

0. 18090268 

*“> 

—0 . 00805 1 1 672 

34 

0.087708324 

T 

—0 . 0056672352 

35 

-0. 0000000596 

4 

-0.0034138164 

36 

-0. 0000000298 

5 

-0.0051245769 

37 

-0.0062124208 

6 

-0.0076082028 

38 

0.01463151 

7 

-0. 010241373 

39 

0.012164063 

8 

-0. 0087890709 

40 

O. 010078562 

9 

-0. 0060926247 

41 

-0. 0068010194 

10 

—0. 0029096406 

42 

-0.01 1 16678 

1 1 

—0. 0039897766 

43 

-0.013444264 

12 

—0. 0069566377 

44 

-0. 0057185646 

13 

-0. 010872209 

45 

0 . 000825 1 369 

14 

-0. 01031487 

46 

0. 0013090034 

15 

-0. 0069459821 

47 

-0. 0055334438 

16 

-0.0019447953 

48 

-0.01199099 

17 

— 0 . 00 151 3800 1 

49 

-0. 014283719 

18 

-0.005201 1646 

50 

-O. 0092434743 

19 

-0.01054764 

51 

-0.0031738635 

20 

-0.0063482042 

52 

O. 0003896314 

21 

-0.0064172242 

53 

-0. 0030636846 


0. 0077392296 

54 

-0. 0086372625 

'“i “T 
JL 

0. 009096365 

55 

-0. 013200269 

24 

0.013419379 

56 

-0.01 1261439 

2j 

-0. 0045658275 

57 

-0.0063766045 

26 

0.0031013507 

58 

-0.001 1074292 

2? 

0. 001939863 

59 

-0.0016479371 

26 

0 . 086Bo 1 067 

60 

-0.0057100602 

29 

0. 17799297 

61 

-0.011455595 

30 

0.94892192 

62 

-0. 01 2233e58 

31 

1.7708108 

63 

—0 . 009 1 46275 

32 

0.951437 

64 

—0 . 003 1 696965 
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TABLE 5. lb 


CHANNEL 

AMF L 1 1 UDE 

CHANNEL 

1 

-0.0007267897 

33 

2 

0 . 00036362 1 3 

34 

*T 

-0.0007288065 

35 

4 

0.000359321 

36 

5 

-0.0007352171 

37 

6 

0.0003502349 

38 

7 

-0.0007465136 

39 

8 

o • 000335336 1 

40 

9 

—0 . 0007638502 

41 

10 

0.0003129798 

42 

1 1 

-0.0007891749 

43 

12 

0.0002799801 

44 

13 

-0.0008258244 

45 

14 

0. 0002303424 

46 

15 

-0.0008794409 

47 

16 

0. 0001508033 

48 

17 

-0. 000961034 1 

49 

18 

— 0 . 0000059596 

50 

19 

-0.0006487506 

51 

20 

0. 0061890148 

52 

21 

-0.0004842821 

53 

22 

0.0089730266 

54 

n 

-0. 0004967653 

55 

24 

0. 01455773 

56 

25 

-0.00047 1 1586 

57 

26 

0.02819664 

58 

27 

-0. 0000016817 

59 

28 

0. 078081675 

60 

29 

0.016213726 

61 

30 

0. 85206997 

62 

31 

1 . 9684365 

63 

T O 

0. 85207033 

64 


AMPLITUDE 
0. 01 61 14548 
0.078082889 
0. 0000000596 
0.028198689 
-0. 0004686741 
0.01 4560707 
-0.0004932102 
0.0089769382 
-0.0004799398 
0 . 006 193867 
-0.0006433618 
O. 0000000065 
-0.0009544246 
0.0001579961 
-0.0008715394 
O. 0002389747 
0. 0008 163354 
O. 0002902792 
-0. 0007779046 
0.0003252768 
0.0007504617 
0.0003499581 
0.0007305867 
O. 000367618 
0.0007162648 
O. 0003801 139 
-O. 0007061326 
0.0003887031 
-O . 000699407 4 
0.0003940842 
-0 . 0006955029 
0.0003967483 
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TABLE 5. 16 


CHANNEL 

AMPLITUDE 

1 

-0. 00077 1 7316 

2 

0. 0003937155 

“7 

“0. 0007738003 

4 

0. 0003893094 

c~ 

wJ 

-0.0007804257 

6 

0. 0003798906 

7 

-0. 0007920204 

e 

0.0003645715 

9 

-0.0008099162 

10 

0.00034 15833 

1 1 

-0. 0008359729 

1 2 

0. 0003076236 

13 

-0.0008735884 

14 

0.0002566777 

15 

-0.0009285968 

16 

0.0001751472 

17 

-0. 0010121401 

18 

0. 000014767 1 

19 

-0.0006910345 

20 

0. 0063565653 

21 

— 0 . 0005224662 

'"'i 

0. 0091902679 

23 

—0. 0005352072 

24 

0. 014859481 

OC 

/ v l 

-0.00051 19497 

26 

0. 026640091 

r> “7 

jL. t 

-0. 00004097 1 6 

26 

0. 07860754 

29 

0.016175929 

30 

0.85347605 

31 

1 . 9684162 

“TO 

0. 85596037 


CHANNEL 

AMPLITUDE 

TIT 

0.017002968 

34 

0. 12062502 

35 

0. 0983834 So 

36 

0.070736912 

37 

0. 0003068298 

38 

0.016410794 

39 

-0. 0005344908 

40 

0.010318644 

41 

-0.0005498171 

42 

0. 0063389936 

43 

-0. 0007362581 

44 

-0. 0000175606 

45 

-0.0010287495 

46 

0.0001623952 

47 

-0. 0009356655 

48 

0.000252703 1 

49 

-0. 0006743219 

50 

0. 0003093 1 

51 

-0. 0000320084 

52 

0. 0003476706 

53 

-0. 000801898 

54 

0. 0003746022 

55 

-0. 0007801768 

56 

0. 000393806 

57 

-0.0007646617 

58 

0.000407364 

59 

-0. 000/53648 

60 

0.0004166673 

61 

-0. 0007463681 

62 

0.0004224689 

63 

-0.0007421259 

64 

0. 0004253136 
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TABLE 5.17 


CHANNEL 

AMPLITUDE 

1 

-0. 0007692082 

'—i 

0. 000400869 

3 

-0. 000771243 

4 

0.0003965298 

5 

-0. 0007777035 

6 

0.0003873242 

/ 

-0. 0007890451 

6 

0. 0003724316 

9 

— 0. 0008064367 

10 

0 . 00035002 

1 1 

-0. 0008317996 

12 

0 . 0003 1 6902 

13 

-0. 000868535S 

14 

0. 0002672436 

15 

-0. 0009220879 

16 

0.0001877573 

17 

-0.0010035612 

16 

0. 0000314856 

19 

-0. 000684932 

20 

0 . 0063850 1 <1*7 

21 

-0.0005197326 

'“*> 
a- *1 

0.0092325676 

*7 

-0. 0005312714 

24 

0. 014927295 

2j 

-0 . 0005066058 

26 

0. 02876861 

27 

-0. 0000407096 

28 

0.079146251 

2S 

0. 015836433 

30 

0. 85279739 

31 

1.968987 


0.85529232 


CHANNEL 

AMF'L I TUDE 

-r-r 

0. 016646266 

34 

0. 12092602 

35 

0. 098412037 

36 

0.070831031 

37 

0. 0002950793 

38 

0. 018493826 

39 

-0.0005306772 

40 

0. 010366028 

41 

-0.0005469068 

42 

0.0068696984 

43 

-0.0007297598 

44 

-0. 0000000289 

45 

-0.0010197775 

46 

0. 0001753429 

47 

-0. 0009290504 

48 

0. 0002633152 

49 

-0. 0008692722 

50 

0. 0003185329 

51 

-0.00082797 

52 

0. 0003559383 

53 

-0.0007986321 

54 

0.0003822021 

55 

-0.0007774544 

56 

0. 000400921 1 

57 

-0. 0007623264 

58 

0.0004141744 

59 

-0.0007515688 

60 

0.0004232136 

61 

-0.0007444961 

62 

0.00042887 1 

63 

-0. 0007403672 

64 

0. 0004316319 
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TABLE 5. 16 


CHANNEL 

AMPLITUDE 

CHANNEL 

1 

—0.001023 5 4 5 9 

77 


0.0001517312 

34 

“7 

-0 . 00087 1 5268 

35 

4 

0 . 0003 1 25338 

36 

5 

-0.0008565909 

37 

6 

0. 000 1 587073 

38 

7 

-0. 0010444 1 12 

39 

8 

0. 0001074481 

40 

9 

— 0. 0009093343 

41 

10 

0. 00030038 

42 

1 1 

-0. 0008516007 

43 

12 

0 . 000 1 363355 

44 

13 

-0.001 1182111 

45 

14 

-0 . 0000 1 39704 

46 

15 

—0 . 00 10147353 

47 

16 

0. 0002301 898 

48 

17 

-0. 0008770221 

49 

18 

0.0000143998 

50 

19 

-0. 000923537 

51 

20 

0. 0060222498 

52 

21 

-0.0006596057 

53 

Hi 

0. 0094446745 

54 

23 

-0. 0001 106963 

55 

24 

0. 015344099 

56 

25 

-0. 0005510989 

57 

26 

0. 028871795 

58 

27 

0. 0000672724 

59 

28 

0 • 0823620 1 9 

60 

2V 

0.017755866 

61 

30 

0.84716296 

62 

31 

1 . 9684452 

63 

32 

0.84725189 

64 


AMPLITUDE 
0. 017859034 
0. 082392663 
0 

O. 02876374 1 
-0. 0006075483 
0.015388817 
-O. 0000001227 
O. 009529734 
-O. 0006703536 
0.0059259832 
-0.001022293 
-0. 0000000149 
-0. 0007902887 
O. 0003495659 
-0.0009598439 
-O. 000067477 
-0.0012323714 
0.000062424 
— O. 0008 1 54805 
0.0004247015 
-0. 000793864 
0.0001222414 
-O. OOl 138442 
0. 00004 1 1 073 
-O. 0008865158 
0.0004072865 
-O. 0007151581 
O. 0002504764 
-0. 0010566916 
0. 0000239295 
O. 0009726973 
0. 0003452739 
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TABLE 5. 19 


CHANNEL 

AMPLITUDE 

1 

0. 0 1 4027503 

c 

-0. 014028227 


0. 014034309 

4 

-0.014041785 

5 

0. 01 40551 65 

6 

— 0 • 014 0 7 03 3 5 

/ 

0. 014092386 

6 

-0.014116675 

9 

0.014149866 

10 

-0.01418641 

1 1 

0.014234834 

12 

-0. 014288335 

13 

0. 014359378 

14 

-0.014438513 

15 

0.014545067 

16 

-0. 014665774 

17 

0. 014832503 

18 

-0. 01502596 

19 

0. 015304765 

20 

-0.015639521 

21 

0.01615287 

22 

-0. 016800832 


0.017695328 

24 

-0.019385239 

25 

0. 022352934 

26 

-0.026902139 

27 

0. 0397986 1 V 

28 

-0.057549331 

29 

0. 23478037 

30 

1 . 597302 

31 

0.2275579 

32 

-0. 036255151 


CHANNEL 

AMFLITUDE 

C'C' 

—0 . 0000000 121 

34 

0. 190361o 

35 

1.6371415 

36 

0. 19035983 

37 

0. 0000033304 

38 

-0. 03625989 

39 

O. 22756466 

40 

1.5972935 

41 

0.23479062 

42 

-0.057561494 

43 

0.039812744 

44 

-0.026918059 

45 

0 . 02237 1 1 1 

46 

-0. 019405698 

47 

0.017918076 

48 

-0.016825864 

49 

0. 016180804 

50 

-0. 01567018 

51 

0 . 0 1 5338654 

52 

-0. 0150631 19 

53 

0.014873258 

54 

— 0 .01471 0266 

55 

0. 014594066 

56 

-0. 014492093 

57 

0. 014418321 

58 

-0.014352719 

59 

0.01430596 

60 

-0.014264372 

61 

O. 014236153 

62 

-0.01421 174 

63 

0.014198057 

64 

-0.014187049 



TABLE 


5 . 20 


CHANNEL 

AMPLITUDE 

1 

0. 013326434 

o 

A. 

-0. 013327075 

~z t 

0. 013332912 

4 

-0 . 0 1 334006 

5 

0. 013352912 

6 

—0. 013367363 

7 

0.013388529 

8 

-0.013411701 

9 

0.013443629 

10 

-0.013478592 

1 1 

0.01352516 

12 

—0. 013576454 

13 

0. 0 1 3644673 

14 

—0. 013720875 

15 

0. 013823727 

16 

-0. 013939913 

17 

0. 014101386 

18 

—0. 01426B158 

19 

0. 014559092 

20 

-0.01486351 1 

21 

0. 015384625 

22 

-0.016015731 

23 

0.017089425 

24 

-0.018549023 

25 

0.021477759 

26 

-0. 025964499 

27 

0.038774937 

28 

—0. 056363199 

29 

O. 23334444 

30 

1 . 5994115 

31 

0.22451079 

32 

-0.024277404 


CHANNEL 

AMPLITUDE 

33 

0. 07948482 

34 

0. 20233938 

35 

1 . 6340945 

36 

0. 19246924 

37 

-0.0014325902 

38 

-0. 035073802 

39 

0.226541 

40 

1.5982311 

41 

0.23391546 

42 

—0. 056725346 

43 

0. 039006829 

44 

-0.026132951 

45 

0.021602863 

46 

—0.018649761 

47 

0.017172426 

43 

-0. 01608807 

49 

0. 01544968e 

50 

-0.01494433 

51 

0. 014617335 

52 

-0.014345434 

53 

0.014158753 

54 

—0.013998368 

55 

0.013884398 

56 

-0.013784283 

57 

0.013712084 

58 

-0.013647805 

59 

0.013602111 

60 

—0. 01356 1 386 

61 

0. 013533884 

62 

-0.013510015 

63 

0.013496719 

64 

-0.013485905 
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APPENDIX 1 


DIMENSION A (20), Q(20), W(257) , F(513), FT (513), X (257) 
PRINT 10 

10 FORMAT ( ENTER NO. OF PARTITIONS ') 

READ * , N 
PRINT 20 

20 FORMAT ( ENTER NF ') 

READ * , ITT 
PI =3. 141592/ 

T —0 . 

PRINT 30 

30 FORMAT ( ENTER NO. OF WINDOWS ' ) 

READ * , M 
DO 50 1 — J . I I 

PRINT * , 'WINDOW NUMBER ,1 
PRINT 43 

43 FORMAT < ' ENTER CuEF . ) 

REAij * , A ( I ) 

PRINT 44 

44 FuKHhT ( ENTER SUBSCRIPTS > 

0 READ * , Q ( I ) 

1 NF— 2 

2 NF — NP*2 
PRINT 'NR- , Nr 
IF(NF.LT.N) GO TO 52 
NF2-2+NF' 

DO 99 K— 1 , NF 
WOO — 0. 

IF(K.GT.N) GU TO 99 
DO 90 1 1 - 1 , M 
I F ( T . EQ . 0 . ) W ( K ) — 1 . 

IF (T . EQ. u. ) GO TO 90 
I F ( O < I I ) . G ( . 1 . ) GO T 0 7 0 
IF < T . LE . NF ) W ( K ) =W ( K ) + A ( I I > 

70 W2- ( 1 . -4. * (T**2) ) 

IF <0 ( 1 1 ) . EQ. 2. ) W GO =WOO +A (II) *W2 
I F < Q ( 1 1 ) . EQ . 3 . ) W ( K ) = W ( K ) + A (II) »COS ( P I *T ) 

IF (Q ( 1 1 ) . EQ. 4. ) WOO=WOO +A (II)* (SIN (2. *PI*T) / (2.*PI*T) ) 
I F < Q ( 1 1 ) . EQ . 5 . ) W ( K ; W ( K ) +A ( I I ) * ( 1 . +2. * ( ABS (T) ) ) 

IF (Q ( I I ) . EQ. 6. ) W (K> — W <K ) + A ( 1 1 ) * ( . 5+. 5*C0S (2. *F I*T ) ) 

Wl-1 . 

IF (Q ( 1 1 > . EQ. 7 . ) W GO -W GO + A ( I 1 ) * ( . 54+ . 46*C0S (2. *PI*T) ) 

IF <Q ( 1 1 ) . EQ. 1 . ) W GO =W GO +A ( 1 1 ) *W1 

W8= ( 1 . /PI ) *ABS (SIN (2. *F'I*T ) ) + ( 1 . -2. *ABS <T ) *C0S<2. *PI*T ) ) 
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IF <Q < I I > . EG. B. > W <F)=W(K) +A (II) *Wb 
W9-. 42+ . 5*C0b < 2. *P J *7 ) + . Ofe*COS ( 4 . *Y 1*7 » 

IF<Q< I I ) .EG!. 9. ) W (K)=W (K> +A(II) +W 9 

IF (T. 6 1 . < .25) ) GO Id 71 

IF (Q ( I I ) .61 . 10. ) GO TO B9 

W10- 1 . - <24. *ABS <7 ) **2 ) * < 1 . -2. + AB S < T ) ) 

IF(0(IU . EG - 1 0 > W C K ) =W ( K ) + W 1 0 
71 I F < G* ( 1 I ) . LT . 1 O . i GO TO 90 

I F ; ( Q ( II). GT . 1 0 . ) GO TO 89 

I F < 1 . GE . \ . 25 > ) W v K ) - W ( K ) -+-A (II) *2* ( 1 . -2 . * ABB ( T ) ) **3 
ft 9 IF '>0(11;. EQ. 11.) W<K> < K ) +A ( 1 1 ) * ( 1 . -2. * (ABB (T) ) > 
IF (1 . GE. < . 15) ) Gu TO Sfe 
iX ( Y ) (K> + A ( I J. ) *Ui 1 

GO ‘I Ci 90 

68 W 1 2= ; < . 5+ . 5+ CDS y ( 2*F‘ I * ( ABB v T ) - . 15) ) / . 7 ; ) 

I F ( 0 v I I ) . EQ . 12 .) U) ( K ) ~W < K ) + A ( 1 1 ) *W 1 2 

90 CONTINUE 

IF (K. EQ. 1 > GO TO 99 
W (NF'2-K+2) =W U.) 

9 9 T — T + 1 . / ( T h * 2. > 

DU 80 L==1,NF2 
W ( NP+ 1 ) = 0 . 

FT ( 2*1. ) =-0 . 

80 r T < 1>. L- 1 t = W ( L ) / ( TM*I * 

DO 101 K2“ 1 , 2-*TJP 2 
101 X<K2)— 0. 

X (37 ) = . 000O975 
X ( 39 ) — • 0033 475 
X (4 1 ) — . 0000965 
X ( 43 ) - . 0050 1 95 
X ( 45 ; — . 0000955 
X ( 47 > == . 0083305 
X (49) 0000945 
X < 5 1 ) — . 0163905 
X ( 53 ) = . 0000935 
X v 55 ) =. 0456593 
X < 57 ) - . 00009 i 3 
> (59) ^.41 15925 
X ( 6 1 ) = 1 . 0 
X (63) -X (59) 

X (65) =X (57) 

X (67 ) X ( 55 ) 

X <69) =X (53) 

X (71 ) =X <51 ) 

X (73) =X (49) 



X <75; --X (47) 

X (77) =X (45) 

X (79 ) -X (43) 

X (81 ) = X (41 ) 

X (83) =X (39; 

X (85) =-X (37) 

X ( 22 1 / -X (37 ) 

X (219) -X (39) 

X ( 2 1 7 ) = X (41) 

X < 2 1 5 ) — X (43) 

X (213) =X (45) 

X (21 1 ) =X (47) 

X ( 209 ) - X ( 49 ) 

X <207 ; = X (51 ) 

X (205) — X (53) 

X ( 203 ) =X ( 55 ) 

X (201 )=X (57) 

X < 199)=X (59) 

X ( 1 97’ ) = X (61) 

X ( 195)=X (63) 

X ( 193) =X (65) 

X ( 191 ) =X (67) 

X ^ 1 Q9 > -X ( 69 ) 

X * 3 87 ) =X (7 1 ) 

X ( 1 85 ) —X (73) 

X * 1 83) -X (75; 

X ( 1 8 1 ) =X ( 77 ) 

X ( 179) -X (79) 

X U77)=X (81 ) 

X (175) =X (83; 

X ( 1 7 3 ) = X ( 85 ; 

502 CQN7 3 NUfc 

CALL FFT (NF 2 , 1 , X > 

DO 102 K3= 1 . NP2 
F r <2*K3> =0 . 

102 FT (2*K3-1 ) =W <K3> *X (2*L3-1 ) 

WRITE (75,405) (k , w (k) , k=l ,np2) 

405 FORMAT ( 14 ,2X ,F16. 8) 

WRI TE (76,*; (FT (K) ,k>l ,2*NP2) 
CALL FFT (NF'2,-1 ,FT) 

DO 401 IK= 1 , 2*NF“2 

401 FT ( IK) =FT < IK) /t>4. 

205 FORMAT <F4. 1 ,2X,Ele>.8> 

DO 206 I A- 1 , 2*NF‘2- 1 , 2 

206 X ( IA) = ( IA+1 ) /2 



■'J Q rj] 4^ cj r o 


146 

201 


WR 1 7 E < 5o , 205 ) <X(IA),FT(IA),I A— 1 
WR I TE ( 5 1 , * ) ( FT ( I A) , I A-2 , 2+NP2 , ? ) 

WRITE (59 ,*> < < 1 A+ 1 ) / 2 , I A-65 , 77 0 ) 

PRINT 146 

FORMA I ( ' REAL IN 50, IM IN 51 ' i 

STOP 

END 


127,2) 
(FT < IA) 


1 A-65 , 77 , 


SUBROUI INE FF-7 (NN, I SIGN, DATA) 
DIMENSION DATA (600) 

N=2*NN 

J=1 

DO 5 1=1, N, 2 
IF(I-J) 1 ,2,2 
1 TEHF'R- DA TA ( J > 

TEMP1=DATA ( J+l ) 

DATA ( J ) =DATA ( I ) 

DATA (J+l )=DATA< 1+1) 

DATA ( 1 ) =TEMFR 
DATA ( 1+1 ) =TEMPI 

M-N/2 

IT ( J-li) 5 ,5,4 
J = J — M 

M-M/2 

IF (M— 2) 5,3,3 

J —iJ +H 

MM A »2 

I F (MMAX-N) 7,10,1 0 
I BTEF‘-2*MMAX 

THETa-6. 263 1853 /FLOAT ( ISIGN*MMAX) 
SINTH=SIN (THETA/2. ) 

W B T P R - - 2 . * S 1 NT H * S 1 N T H 
WSTPI-SIN (THETA) 

WR= 1 
W I =0 

DO 9 M— 1 , UMAX, 2 
DO a I=li , N, I STEP 
J = I + MM A X 

TEMPR-WR*DATA ( J ) -W] *DATA < J + l ; 

TEMP I =WR+DATA ( J + l ) +WI*DAT A (J ) 

DAT A ( J ) =DA7 A ( I ) - I EMFP 
DATA ( J+l ) =DATA ( 1+1 ) -TEMPI 
DAT A ( I ) =DA1 A(I) +TEMPR 
DAT A ( I + 1 ) =DAT A ( I + 1 ) +TEMF 1 
TEMPR=WR 

WR=WR*WBTPR-WI*WSTPI+WR 
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v W I W 1 *WSTPR+ T h.MF'R*WSTF' I +W 1 

MMAX=ISTtp 
(30 TO 6 

10 RETURN 

LND 
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